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Abstract 

Properties of the pure solitonic r-function and potential of the heat equation are 
studied in detail. We describe the asymptotic behavior of the potential and identify 
the ray structure of this asymptotic behavior on the a;-plane in dependence on the 
parameters of the potential. 



1 Introduction 

The Kadomtsev-Petviashvili (KP) equation was derived as a model for small-amplitude, 
long- wavelength, weakly two-dimensional waves in a weakly dispersive medium [Tj. This 
equation is a (2+l)-dimensional generalization of the celebrated Korteweg-de Vries (KdV) 
equation, and from the beginning of the 1970s |H|3] it is known to be integrable. There 
are two inequivalent versions of the KP equation: KPI and KPII. Here we consider KPII 
equation 

(lit - 6uu Xl + u XlXlXl ) Xl = -3u X2X2 (1.1) 

(the KPI equation has an opposite sign in the r.h.s.), where u = u(x,t), x — (2:1, £2) and 
subscripts x\, xi and t denote partial derivatives. The KPII equation is integrable since 
it can be expressed as compatibility condition [C,T] = of the Lax pair C and T, where 
operator 

C(x,d x ) = -d X2 +d^-u(x) (1.2) 
defines the well known equation of heat conduction, or heat equation for short, and 



T(x,d x ,d t ) = d t + 4d* 1 -6ud x 



3u, 



3d Xl u x 



(1.3) 



The spectral theory of the operator (|1.2I) was developed in [4HZ] m the case of a real 
potential u(x) rapidly decaying at spatial infinity, which, however, is not the most interesting 
case, since the KPII equation was just proposed in OQ in order to deal with two dimensional 
weak transverse perturbation of the one soliton solution of the KdV. In fact, if u 1 (t, x-^) 
obeys KdV, then u(t, x 1 , x 2 ) — Ui(t, x 1 + /ix 2 — 3[i 2 t) solves KPII for an arbitrary constant 
/j£R. In particular, KPII admits a one soliton solution of the form 



m(x, t) 



scch 



-Xi 



-x 2 - 2(k\ - nl)t 



(1.4) 



KPII 



heatop 



Lax 



1-sol 



where ki and K2 are real, arbitrary constants. 
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A spectral theory of the heat operator (|1.2I) that also includes solitons has to be built. In 
the case of a potential u(x) rapidly decaying at spatial infinity, according to [4-7 , the main 
tools in building the spectral theory of the operator (11.21) , as in the one dimensional case, 
are the integral equations defining the Jost solutions. However, if the potential u(x) does 
not decay at spatial infinity, as it is the case when line soliton solutions are considered (see, 
e.g., (jl.4[> ). these integral equations are ill-defined and one needs a more general approach. 
In solving the analogous problem for the nonstationary Schrodinger operator, associated 
to the KPI equation, the extended resolvent approach was introduced [5]. Accordingly, a 
spectral theory of the KPII equation that also includes solitons has to be investigated using 
the resolvent approach. In this framework it was possible to develop the inverse scattering 
transform for a solution describing one soliton on a generic, rapidly decaying background [5], 
and to study the existence of the (extended) resolvent for (some) multisoliton solutions [10] . 

However, the general case of N solitons is still open. In particular this is motivated by a 
complicated asymptotic behavior of the pure solitonic potential on the x-plane. Here we give 
a detailed description of this behavior in terms of the soliton parameters. The general form 
of the multisoliton potential was derived in |llj by means of the rational transformations of 
the scattering data and in [TU] by means of the twisting transformation. In jT^] we proved 
that results of both these constructions coincide one with another and also with the form of 
the multisoliton potential given in terms of the r-function by Biondini, Kodama et al, see 
review article |13| and references therein. 

The paper is organized as follows. In Sec. 2 we write down the multisoliton potential. 
They are labeled by the two numbers (topological charges), N a and Nf,, that obey condition 

N a ,N b >l, (1.5) |nanb 

and showing at large spaces Nb "ingoing" rays and N a "outgoing" rays. By using notations 
of [12] . we give for these potentials two dual representations in terms of r- functions, which 

depend on 

Af=N a + N b , (1.6) |Nnanb 

(so that Af > 2) real parameters K n and on a iV a x N b constant matrix d. In Sec. 3 we 
study the asymptotic behavior at large x of the multisoliton potential in details and show 
that the round angle at the origin can be divided in M angular sectors such that along the 
directions of their bordering rays the potential has constant soliton like behavior, while along 
directions inside the sectors the potential has an exponential decaying behavior, which we 
derive explicitly. In particular, for a special subclass of the solitonic potentials we present 
explicitly the ray structure in terms of the parameters K n . 

2 Multisoliton potentials 

Let we have Af real parameters 

Kl < K2 < • • ■ < (2.1) kappas 

and introduce the functions 

K n (x) = K n x x + «£x 2) n=l,...,M. (2.2) [kpT 

Let 

e K ^ =diag{e K "^}^ =1 , (2.3) 
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be a diagonal TV xTV-matrix, V a TV XTV& constant real matrix with at least two nonzero 
maximal minors, and let V be an "incomplete Vandermonde matrices," i.e., the Nj, x TV- 
matrix 

/ 1 ... 1 



v = 



K N b -l 



Then, the soliton potential is given by 

u(x) = -2^ 1 logT(a:), 
where the r-function can be expressed as 

t{x) =&et(Ve K{x) V), 



(2.4) 



(2.5) [ux] 



(2: 



tan 



see the review paper [13] . references therein, and [12] where the same notations have been 
used. 

There exists (see [14], [TO], and [12]) a dual representation for the potential in terms of 
the r-function 

(2.7) 



r'(a;)=det (t> r e~ K ^'yV' > \ 



where T> is a constant real N a x TV- matrix that like the matrix V has at least two nonzero 
maximal minors and that is orthogonal to the matrix T> in the sense that 



V'V = 0, 

being the zero in the r.h.s. a N a x ^-matrix, and where V' is the TV xNb- matrix 

1 . . 



N a -1 



v' = 



(2.8) 



(2.9) 



tail' 



d!2 



sol ' 



1 



N a -1 
V 



and 7 the constant, diagonal, real TV x TV-matrix 



AT 



7 = diag{7„}^ =1 



(2.10) 



gamma 



This dual representation follows if one notices that 



t(x) = (_1)^+JW- 



i)/2 (n eK " (x) ) 



K(/si,...,/s JVa+ jvjT / (a;), (2.11) 



tautau 



where V denotes the Vandermonde determinant 



V(fti,...,Kjv) = J][ (K„ - Km). 

l<m<n<J\[ 



(2.12) 



Matrices 2? and V ' obey rather interesting properties [12] . Taking into account that the 
matrix T> has nonzero maximal minors we can always write it in the form 



V = 7T 



d 2 



(2.13) 



blockl 
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where tt is a AT x A/"-permutation matrix such that the Nb x iV^-matrix d 2 is nonsingular. 
Then, we can rewrite this equality in the form 



V = 7T f J | d 2 (2.14) I block 

V h N b J ' 1 

where c? = d^^ 1 is a constant real N a x iV^-matrix and IJ/v,, the Nb x iVb unity matrix. 

Now, we notice (see [13]) that the expression for the potential in (|2.5[) . as well as condition 
(|2.8p . are left unchanged if in (12.61) and (|2.7[) . respectively, we multiply the A/" x iV^-matrix 
2? from the right by a nonsingular constant Nb x A^b-matrix and the N a x TV-matrix T>' 
from the left by a nonsingular constant 7V a x iV a -matrix. We conclude that the matrix d 2 
is unessential and that the matrix T> can be chosen to have a simple block structure. It 
is clear, however, that choice of matrixes tt and d is not unique as the matrix T> can have 
different nonzero maximal minors. 

Let us write now T> tt = (d' 1: d' 2 ), where d[ is N a x iV Q -matrix and d' 2 is N a x iVf,-matrix. 
Then by ()2.8j) we get d! 2 = —d[d, where d is the same N a x Nj, matrix as in (|2.14|) . Thus, 
for V' we get the block structure 

V = d' 1 (E Na ,-d)<K\ (2.15) | block' 

and thanks to equivalence of r and r' we get that det d[ ^ and, then, also V ' can be 
chosen to have a simple block structure. 

Moreover, by means of the block structures given in ()2.14|) and (|2.15|) we get 

2? f V = d\{E Nb + dU)d 2 , V'V' 1 = d' 1 (E Na +dd^)d' 1 \ (2.16) [d!4 

where f denotes the Hermitian conjugation of a matrix (in fact, transposition here). So both 
matrices V and T>' V' are positive and then T> and T>' admit, respectively, left and right 
inverse (see [15]). Precisely, we have 

(V) { - 1) V = E Nbl V'(V')^ 1} =E Na , (2.17) [dlT] 

with 

{V) ( ~ 1] = {V^VY X V\ (V') ( ~ 1} = V'\V'V A Y 1 . (2.18) [dis" 

Products of these matrices in the opposite order give the real Hermitian M x A/"-matrices 

P = V(V) ( ~ 1] = V{V^ V)- 1 (V)\ (2.19) 

p' = (v') { ~ 1} v' = {v')\v'v A )- 1 V , (2.20) [die" 

which are orthogonal projectors, i.e., 

P 2 = P, (P') 2 =P', PP' = = P'P, (2.21) [d22" 
and complementary in the sense that 

P + P' = E M - (2.22) [d23~ 

Orthogonality of the projectors follows from (|2.8[) and the last equality from obvious relations 
of the kind (E Nb + dU)- 1 d^ = d^(E Na +dd^)~ 1 . 
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Finally, let us mention that if all maximal minors of the matrix T> are nonzero, the 
permutation 7r in f|2 . 14[) and (|2. 15|) can be arbitrarily chosen. Then, we can chose it equal 
to the identity, getting the simplified representations 



V = 



d 

E Nb 



da, V = d' 1 {E Na ,-d), 



so that in this special case without loss of generality we choose 



V 



d 

E Nh 



V = {E Na ,-d). 



(2.24) 



In order to study the properties of the potential, it is convenient to use an explicit repre- 
sentation for the determinant (|2.6[) . By using the Binet-Cauchy formula for the determinant 
of a product of matrices we get 



t(x) 



l<n 1 <n 2 <---<n^ 



fn 1 ,...,n Nb 1 



(2.25) 



o <N 1=1 
with coefficients fm,n 2 ,...,n Nb given by 

fm,n 2 ,...,n Nb = V{n ni , . . . , n nNb )V(n 1 , . . . ,n Nb ). 
Here we used notation (|2.12l) for the Vandermonde determinant and notation 

ni ,N b 



V{m, ■ . .,n Nb ) = det 



(2.27) 



for the maximal minors of the matrix T>. Notice that these coefficients are invariant under 
permutations of the indices. 

From (12.25)) it follows directly that condition 



f ni ,...,n Nh > 0, for all 1 < m < n 2 < ■ ■ ■ < n Nb < TV, 



(2.28) 



is sufficient (see [13]) for the regularity of the potential u(x), i.e., for the absence of zeros of 
t{x) on the x-plane. Thanks to (|2.1|) . (|2.12j) and p.26| this condition is equivalent to the 
condition that all maximal minors of the matrix T> are non negative. In |12j it was mentioned 
that condition (|2.28|) is also necessary for the regularity of a potential under evolution with 
respect to an arbitrary number of higher times of the KP hierarchy. In [13] it is suggested 
to decompose soliton solutions of KPII into subclasses, associated to the Schubert cell on 
Grassmanian. It is proved there that condition (|2.28p is necessary for the regularity of the 
all solutions associated to a cell. However, problem of finding of the conditions necessary 
for the regularity of the many soliton solution under evolution with respect to KPII only is 
still opened. 

If the matrix V can be chosen as in (|2.24[) , then introducing new N a x iV^-matrix d as 



d n ,i — dN a +i-n,i(— 1) 
we easily get by ()2.27|) that 

V(m, . . . ,UN b ) = det Hdnj./H, 



■,N a , 1 = 1, 



(2.29) 



(2.23) blockreg 



block2 



tauf 



(2.26) 



Do 



reg 



dd2 
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wherei = !,...,&, 1 = 1,..., Nb, I ^ rij — N a , j = k + 1, . . . , Nf, 



(2.30) 



ddl 



lemmal 



and where number k, 1 < k < Nf, is denned by a condition < N a < rifc+i. This proves 
that if we require condition (|2.28[) on coefficients f ni ,....n Nb , then the matrix d is a totally 
nonegative (positive, if all inequalities in (|2.28l) are strict) one (see ( [T6l[T7] ). 

In what follows it is convenient to continue indexes of n n , f ni ,...,n N > etc, periodically by 
TV on the whole Z by means of condition: 



n n (mod TV). 



Then Eq. (|2.25|) can be written as 
t(x) = 



n<ni<ri2<---<niv. <JV"+n— 1 



n p K„,(x) 
e 



(2.31) 



(2.32) 



where the r.h.s. is independent of n, and where coefficients f ni ,...,n N are defined by (]2 .26[) 
taking (|2.3ip into account. Thanks to (|2.1j) and (|2.3ip we can prove the following lemma. 



Lemma 2.1 Let 



91, 



( K i K m)( K n K n+N b 



K n+N b ) ■ 



Then 



fJi.i 



> /or any n G 



I = n, . . . ,n + Nj), 



b, ■ ■ 



-n, 



(2.34) 



and equality takes place only when I and m independently take values n or n + Nb by modAf. 



Proof. Thanks to (|2.3ip it is enough to prove the lemma for n = 1, . . . ,7V. For these 
values of n the k„'s are ordered as in (|2.ip . For other values the ordering is obtained by 
coming back by means of (|2.3ip to values belonging to the interval n — 1, . . . ,7V. Thus, 
when n = 1, . . . , N a then n, n + Nb and / are not more than TV (see ((US])), but instead of the 
interval for m given in (|2.34[) wc get two intervals: m = n + Nb, ■ ■ . ,M and m — 1, . . . ,n. 
Then for given values of n, I and m in the first interval lemma follows because the first factor 
in (|2.33p is nonpositive, the second one negative and the third one nonnegative. Zero of the 
first factor is possible only for I = m and then I — m — n + Nb and of the third factor only 
when both / = n and m = n+Nb (all equalities by mod A/"). In the second case the first factor 
is nonnegative (zero only at I = m = n (mod TV)), the second one is negative and the third 
one is nonpositive (zero only at I = n + Nb, m = n). In the interval n — N a + 1, . . . , J\f we 
write by (|2.3ip K n +N b = Kn-N a > so that again n and n — N a belong to the interval 1, . . . , TV. 
Now the interval of I given in fl2.34|) decouples in two intervals n, ... ,7V and 1, . . . ,n — N a , 
while m = n — N a , . . . , n, so it is positive and less than TV. Then for these values of n, m 
and I in the first interval thanks to (I2.ip the first factor in (12.33)) is nonnegative (zero only 
when I = m = n (mod TV)), the second factor is positive, and the third one is nonnegative 
(zero only when / = n and m = n + Nb). In the case of the second interval the first factor 
is nonpositive (zero only when I = m = n — N a = n + Nb (mod TV)), the second factor is 
positive and the third factor is nonpositive (zero only when I = n — N a = n + Nb (mod TV) 
and m = n). ■ 



modN 



taufn 



(2.33) |TlTl 



112 



G 




Figure 1: Sector a n corresponds to N a + 2 < re < A/" — 1. fig 

3 Asymptotic behavior of the r-function and potential 

u(x) 

Here we study the asymptotic behavior of the function t(x) that by (|2.25p is determined by 
interrelations between functions K n (x) for different n. Since by (|2.2[) 



Ki{x) - K m (x) = (m - K m )(xi + (kj + K rn )x 2 ), for any I, m G Z, (3.1) | linen 

that is, these differences are linear with respect to the space variables, the asymptotic 
behavior must have sectorial structure on the x-plane. In order to describe this structure 
we introduce rays r n given by 

r„ = {x : xi + (nn+N b + n n )x 2 = 0, (n n +N b ~ K n )x 2 < 0), U=1,...,jV, (3.2) | rayn 

that intersect the lines x 2 = ±1, respectively, at the points 

(K n+Nb + K n , -1) for re = 1, ... , N a 

(-K n+Nb - K„, 1) for re = N a + 1, . . . ,Af. (3.3) yOn 
Thanks to ((5TTJ and (|2"3Tj) . 

K„+Ar b - K„ > 0, n = 1, . . . , N a , (3.4) ineql 

K„+Ar b - Kn < 0, n = N a + 1, . . . ,Af, (3.5) ineq2 

therefore, increasing re from 1 to A^ a the ray rotates anticlockwise in the lower half re- 
plane, crosses the positive part of the xi-axis in coming to re = N a + 1 and, then, rotates 
anticlockwise in the upper half ir-plane up to n = Af and, finally, crossing the negative part 
of the xi-axis, for re = Af+1 comes back to the ray r% thanks to (|2.3ip . see Fig.Q] 

Let us assume that some rays, say, r m and r„, where for definiteness m < n, are parallel. 
By (|331 this means that n m+Nb + n m = n n+Nb + n n , i.e., that n n - K m = n m +N b - K n+Nb , 
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where the l.h.s. is positive thanks to (|2.1[) . Then, because of (|2.31[) it is easy to see that the 
r.h.s. can be positive only if 1 < m < N a < n < TV, so that by (|3.2p ray r m is in the bottom 
halfplane and r n in the upper one. Thus rays can be parallel only if they belong to different 
halfplanes, but in the case N a ^ Nb this is possible for a special choice of the parameters n n 
only. On the contrary, in the case N a = Nb all pairs r n and r n+ N a , n = l,...,N a , and only 
these pairs give parallel rays. In the special case N a = Nb — 1 we get two rays producing 
the straight line x\ + (k\ + k^)x2 = that divides the x-plane in two halfplanes. 
Now we introduce sectors o~ n , which are subsets of the x-plane characterized as 



{x : K n -x(x) < K n+ N b -i(x) and K n (x) > K n+Nb (x)}, n = 1, . . . ,7V. 



(3.6) 



Thanks to (|3. 1[) and the discussion above it follows that the sectors a n are sharp (for TV > 2) 
angular sectors with vertexes at the origin of the coordinates bounded from the right (looking 
from the origin) by the ray r n -i and from the left (looking from the origin) by the ray r n . 
With the increasing of n sectors a n are ordered anticlockwise, starting "from the left" with 
the sector oi, that includes the negative part of the xi-axis, and, then, with the sectors 
cr„ (n = 2, . . . , N a ) in the bottom half-plane, the sector <JN a +i "to the right", that includes 
the positive part of the Xi-axis, and the sectors a n (n = N a + 2, . . . ,7V) on the upper 
half-plane, finishing with the sector a^f tangent to the sector o~\, covering in this way the 
whole x-plane with the exception of the bordering rays r n . Therefore, the sectors a n define 
a TV-fold discretization of the round angle at the origin, with n playing the role of a discrete 
angular variable, see Fig. [TJ It is clear that in the study of asymptotic behavior of r- 
function we consider x-plane as vector space, as finite part of x when x — > oo is irrelevant. 
For determining the directions of rays r n and sectors cr„, we introduce the vectors 



Vn = (« 



2 

n+Nb 



), 



1,...,AA, 



(3.7) 



that enables us to give the following definition. 



defl' Definition 3.1 We say that x —¥ oo along the ray r n if x — > oo and there exists such 
function a(x) — > +oo that x — a{x)y n is bounded. This will be denoted as x oo. 

We say that x — >■ oo in the sector a n , n = 1, . . . ,7V, if x — » oo and there exist such 
functions a(x) — > +oo and j3{x) — > +oo that x — a(x)y n -\ — /3(x)y n is bounded. This will 
be denoted as x oo. 

Notice that for x oo 

K n (x) — K n+ N b (x) is bounded and (K n +N b — K n )x2 — > — oo, (3-8) 

and for x — ^ oo 

K n +N„-i(x) - K n _ 1 (x) -> +oo, K n (x) - K n+Nh {x) ->■ +oo, (3.9) 

as follows directly from the definition. In fact we can prove a more general statement. 

Iemma3 Lemma 3.1 Let N a ,Nb > 1, and n € Z arbitrary. Then 

1. if x — ^ oo then Ki{x) — K m (x) — > +oo or bounded for any I — n, . . . ,n + 7V& and 
m = n + Nb, . ■ ■ ,7V ' +n, where boundedness takes place if and only if (I, m) — (n, n), 
(n, n + N b ), (n + N b , n), (n + N b ,n + N b ); 



2. if x oo then Ki{x) — K m (x) — > +oo for any I 
n + N b ,...,Af+n- 1; 



. , n + Nb — 1 and m 



sigman 



pointn' 



def 12 



def 11 
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where summation in indexes is always understood by modJV. 



Proof. Thanks to (13. ip and (|3.7p we have that Ki(y n ) — K m (y n ) = gi, m ,n- Then the first 
statement of the lemma follows directly from Definition 13.11 and Lemma 12.11 In the same 
way we get that Ki(x) - K m (x) = agi, m , n -\ + /%, m ,„ under substitution x = ay n ^ 1 + j3y n . 
Then the second statement of the lemma follows from the second statement of the Lemma 
I2.1l while intervals for I and m appear as intersections of the corresponding intervals in (I2.34[) 
for n and n — > n — 1. I 

Now we can prove the following Theorem. 



| thl | Theorem 3.1 The asymptotic behavior of t{x) for x — > oo is given by 



'n+JVb-l 



oo : t{x) = (z n + z n+1 e K ^+^- K ^ + o(l)) cxp £ K 3 {x) , (3.10) 



3:232 



(cf. (2.26\) ). Terms o(l) are decaying exponentially. 
Proof. Using representation (|2.32[) we can write 

/ n+N b -l \ 

cxp ( ~ S Ki ( x n T ( x ) = 



N,, 



n+N b -l 

fm,...,n Ni exp I ^2,K nj {x)- ^ K ^ x ) 

n<ni<n2<-- <nN, <A r +?i— 1 V .7=1 l—n 



(3.13) 



(3.11) tauasympt 



/n+N b -l \ 
x oo : r(x) = (z n + o(l)) exp I ^ Ki(x) J , 

for any n € Z, where coefficients z n are defined as 

Zn = fn,n+l,...,n+Nt-l = • • • . K n+N b -l) ' D { n , .-.,n + N b -l), (3.12) QE] 



Let us consider first the asymptotics when x — ^> oo. Then each term K n Ax) either cancels 
with some term in the last sum, when the index rij equals by mod Af one index I in the interval 
{n, . . . , n + Nb — 1}, i.e., an interval of the index I in Lemma 13.11 or belong to the interval 
{n + Nb, . . . ,Af+n — 1}. Thus by Lemma 13.11 all these exponents at infinity are negative 
with the only exception of the case where {n\, . . . , nN b } = {n, . . . , n + Nb — l}(modAT). This 
proves (|3.11[) . Proof of (13.10[) goes in the same way with the only difference that we use the 
first statement of the Lemma 13.11 Statement on the asymptotic behavior of the terms of 
o(l) kind is obvious by construction. ■ 

As we see from p.lOj) and (I3.11[) the leading asymptotic behavior is fixed by this theorem 
only if all 

z n >0, n = l 1 ...,M. (3.14) 

This condition is a sufficient for the regularity of the potential at large x, while specific 
examples show that it is not sufficient for the regularity of the potential in all the x-plane. 
We also see that the asymptotics along the ray direction r n is given by the sum of the leading 
terms obtained for x 00 and x 00. Since the exponential factor cancels out when 
expansion (|3.10p of t(x) is inserted in (|2.5I) . the factor in parenthesis gives the ray behavior 
of the potential u(x) at infinity. Thus, if condition p.!4p is imposed potential u(x) has 



znO 
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exactly Af rays r n , obeying nontrivial asymptotic behavior of onesoliton potential (cf. (|1.4l0 
with parameters K n and n n +N b - Explicitly, taking into account (|2.31|) we get the following 
ray structure and asymptotic behavior of the potential along rays: 

X2 — > — oo, x\ + {n n + K n +N b )x2 being bounded: 

u(x) = -2d 2 Xi \og(z n + z n+1 e K ^^- K ^), n=l,...,N a . (3.15) 

and 

X2 — > +00, x\ + (n n + K n+ N a )x2 being bounded: 

u(x) = -2dl \og(z n+Na + z n+Na+1 e K '^- K ^°^), n=l,...,N b , (3.16) 

i.e., N a asymptotic rays in the bottom half plane and Nf, asymptotic rays in the upper half- 
plane. Asymptotic behavior inside sectors <r n is given by (|3.11l) where the only n-dependent 
term is the exponential factor. Taking into account its linear dependence on x and (|3. 14|) 
we get that u(x) decays exponentially in all directions inside all sectors. 

4 Concluding remarks 

Asymptotic behavior of the r-function and potential derived here is based on condition A3. 14|) 
allowing vanishing or negative values of coefficients in (|2.25p different from z n (see l|3.12|) ). 
Thanks to (I3.14|) we were able to find explicitly that the ray directions are given by the 
sums K n + n n+ N b , with n = 1, . . . , J\f, and to derive the exact asymptotic behavior of the 
r-function on the x-plane. Let us mention that analogous result was derived in [18 for a 
very special subclass of solitonic potentials, given by matrix T> in (|2.6[) equal to a product of 
a diagonal matrix by the matrix transpose to V. As we mentioned above, condition ()3.14p 
does not guarantee (in contrast to (|2.28|) ) absence of singularities of the potential u(x). 
But (|3.10|) and (|3.1ip prove that thanks to (|3.14j) singularities of the potential (|2.5j) cannot 
appear in asymptotics, as inside sectors u(x) — > and along rays r n it is given by p,15[) and 
(|3.16l) . that is also finite thanks to (|3.14[) . Let us emphasize that parameters iV a and Nf, 
giving numbers of rays of the r-function and potential in the bottom and upper half-planes, 
correspondingly, appeared in [10] as numbers of poles of the Jost and dual Jost solutions. 
And, finally, we mention that Theorem l3.1l is also valid when all inequalities in the regularity 
requirement (|2.28[) are strict, i.e., all maximal minors of the matrix T> are positive, as this 
condition is stronger than p.!4p that was enough to prove the theorem. 
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